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Determining the escape fraction of ionizing photons during 
reionization with the GRB derived star- format ion rate 
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ABSTRACT 

The fraction of ionizing photons that escape their host galaxies and so are able to 
ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses 
of the reionization era and early galaxy formation. Studies of the reionization history 
normally suffer from a degeneracy between the unknown values for the efficiency with 
which high redshift galaxies turn mass into stars and the escape fraction of ionizing 
photons. Recent gamma-ray burst (GRB) measurements of the star formation rate 
density during reionization provide the first opportunity to break this degeneracy. 
We confront a semi-analytic model for reionization with the GRB-derived star for- 
mation rate, as well as observations of the Lya forest and the CMB. Assuming that 
UV photons produced in star-forming galaxies dominate the reionization process, we 
show that the escape fraction of ionizing photons from high redshift galaxies is 5% 
[log /esc = —1-35 ± 0.15 (68%)] for our fiducial model. This value is reasonably sta- 
ble against uncertainties in the modeling, including the implementation of radiative 
feedback, the possibility of an evolving escape fraction, and the unknown shape of the 
IMF, which in sum contribute ~ 0.2 dex of additional systematic uncertainty on the 
value of escape fraction. 

Key words: cosmology: diffuse radiation, large scale structure, theory ~ galaxies: 
high redshift, inter-galactic medium 



1 INTRODUCTION 

Star-bursting galaxies and quasars have been the leading 
candidates for the sources of the UV radiation required 
to reionize t he hydrogen gas in th e inter-galactic medium 
(IGM) (e.g., lBarkana fc LoebllioOll ) . The quasar population 
is observed to decline quickly at z > 2.5 and so it is believed 
that it was galaxies that contributed the bulk o f UV pho- 
tons that drove reionization (|Madau et al.iil999|: Fan et al.l 



l2002l : ISrbinovskv fc Wvithj2007l : lBolton fc Haehneltll200'iD . 

However, it has been difficult to make a positive case for 
galaxies as opposed to favoring them by disfavoring the case 
for quasars. The contribution of galaxies to the UV radia- 
tion field is dependent on the star formation rate and initial 
stellar mass function, but is also limited by the fraction of 



ionizing photons that escape their host galaxies. If the es- 
cape fraction is small, then star formation had to be very 
efficient at high redshift in order to reionize the Universe. 
The escape fraction is therefore a critical parameter in stud- 
ies of the connection between high redshift galaxy formation 
and reionization. 

Attempts to determine the escape fraction have been 
dominated by direct observations of relatively low redshift 
galaxies, and by numerical simulation. Direct measurements 
of the escape fraction (/esc) are compficated by the fact that 
the intrinsic number of ionizing photons produced by a spe- 
cific galaxy is unknown. A commonly adopted strategy to 
circumvent this limitation is to compare the fiux observed 
at the Lyman limit to the observed fiu x at a frequency where 
the intrinsic em i ssivit y can be inferred (jLeitherer et al ]|l995l: 
Heckman et al.' 2001; D eharveng et all I2OOII : ISteidel et al.l 
2001; Si ana et al.. .2007'l. The escape fraction at the Lyman 
limit can be then be derived using a model description of the 
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star f ormation history (|Leitherer et alj 1 19951 : ISteidel et alj 
I2001I : [s lana et al.ll2007l ). However the results are currently 
uncertain. At redshifts z ~ 1 — 3, observations have sug- 
gested a br oader range of values for /esc, from a few percent 
to ^ 20% (ISteidel et allbOOll: iFernandez-Soto et al.1 12003 



.. 7,' T V i]l2003l : 

IShaplev et al.ll2006l : ISiana et al.ll2007D . Ilnoue et al.l (I2OO6I ) 

have examined the evolution of the escape fraction in the 
redshift range 2 = — 6 using both direct observations of 
the escape fraction and values that they derive from mea- 
surements of the ionizing background. They find that the es- 
cape fraction evolves from /es c ~ 1 — 10%, i ncreasing towards 
high redshift. More recently iFvnbo et al.l 12009 ) have used 
spectroscopic observations of individual GRB afterglows at 
z > 2 to place a 95% confidence upper limit of 7.5% on the 
escape fraction for ionizing photons on the sightlines of these 
GRBs. 

Theoretical modeling has concentrated on the absorp- 
tion of ionizing photons as they propagate through the inter- 
stellar medium towards the IG M. Here also th e results are 
not yet conclusive. For example. I Wood fc Loebl ^2000,) found 
/esc 1% at 2 ~ 10 suggesting that the escape fraction de- 
creases towards high redshift due to the increased density of 
galact ic disks. More recently, iRazoumov fc Sommer-LarserJ 
l|200m used galaxy formation simulations incorporating 
high-resolution 3-D radiative transfer to show that the es- 
cape fraction evolves from fosc ~ 1 — 2% at z = 2.39 to 
. fesc ^ 6 — 10% at z — 3.6. In agreem e nt wi th iFuiita et al.l 
l|2003h . E azoumov fc Sommer-LarsenI l|2006l ) find that in- 
creased supernova feedback at higher redshift expels gas 
from the vicinity of star-bursting regions, creating tunnels 
in the galaxy through which ionizing photons can escape 
into the IGM. They also find that star formation may occur 
at slightly lower densities at higher redshift, which further 
contributes to the increased /esc towards high er redshift. 

Recently, detailed numerical simulations (jGnedin et al.l 

I2OO7I ) have been presented that predict a value for fesc be- 
tween 1 and 3%, for halos of mass M > 5 x 10^° Mq, 
over the redshift range 3 < 2 < 9. This very low effi- 
ciency of reionization would have profound effects on the 
reionization history. In addition to a small escape frac- 
tion in massive galaxies, iGnedin et al.l (|20o3) further pre- 
dict that halos with Af ^ 5 x 10^" have an escape 
fraction t hat is negligibly small . We point the interested 
re ader to Srbinovskv fc Wvithd l|2008l ) and to the review 
bv lCiardi fc Ferraral l|2005h for a more detailed discussion of 
prior work on estimating the escape fraction. 

A less direct avenue to determining the escape frac- 
tion is through modeling of the ionizing background at the 
end of reionization. Models of the reionization of hydrogen 
and the subsequent Lya forest are subject to a degener- 
acy between the star formation efficiency of gas accreted 
into galaxies and the fraction of ionizing photons that es- 
cape the dense gaseous environment of the galaxy to ionize 
the IGM. In order to break this degeneracy, the true cos- 
mic star formation rate must be known. For example, the 
escape fraction has been estimated at high redshifts by com- 
bining the ionization rate in the IGM with an estimate of 
the star formation rate obtained from measurements of the 
galaxy luminosity function ( e.g., Bolton fc Haehnelt 20071 : 



20081 ). 



ISrbinovskv fc Wvithd l2008l : iFaucher-Giguere et all 
Assuming that the mean- free-path is accurately estimated, 
the escape fraction follows from comparison of the total 



number of photons produced with the number that are ion- 
izing the IGM. The primary drawback of this technique is 
that the observed luminosity functions are subject to uncer- 
tain corrections for extinction in galaxies, as well as a flux 
limit. The latter fact means that an extrapolation to lower 
galaxy luminosities must be performed in order to obtain a 
total star formation rate. 

GRBs are now beginning to probe the star forma- 
tion rate density in the era during which the Universe is 
thought to have b ecome reionized (jSalvaterra et al.] l2009l : 
iTanvir et al.ll2009l ). In this paper we employ t he recent de- 
term ination of the star formation rate density (iKistler et al.l 
I2OO9I ) out to 2 ~ 8.5 based on a compilation of high redshift 
GRBs. As described below, this measurement of the cosmic 
star formation rate with GRBs is not subject to the usual 
large corrections required for extinction and unseen faint 
galaxies. By modeling the reionization of hydrogen we deter- 
mine the star formation efficiency and escape fraction that 
is required to simultaneously reproduce the optical depth to 
electron scattering of CMB photons, the star formation rate 
density during the reionization era, and the ionizing back- 
ground at the end of reionization. Inclusion of the star for- 
mation rate density as a constraint allows the star formation 
efficiency and escape fraction to be separ ately determined. 

O ur an alysis is similar to that of iChoudhurv fc Ferraral 

l|2005l . l2006l ). who modeled the reionization history and com- 
pared model predictions to a range of observables, including 
the star formation rate, the ionization rate from the Lya 
forest and the optical depth to electron scattering of CMB 
photons, as well as the evolution of Lyman limit systems 
and the temperature of the IGM among others. However, 
here we limit our attention to the first three of these ob- 
servables that most directly relate to determination of the 
escape fraction, and which can be most reliably predicted 
by an analytic model. Most importantly, we make use of 
the recent determination of star formation rate density from 
GRBs, which for the first time provide a measurement of the 
total star formation deep into the reionization era. We also 
quantify the statistical and systematic uncertainties on the 
model parameters. 

The paper is organized as follows. In § [2] we discuss 
the determination of star formation rate density from high 
redshift GRBs, and review a simple estimate of the ioniz- 
ing photon budget for reionization. We then describe other 
observational constraints (§|3]), as well as our semi-analytic 
model for reionization (§|31). The resulting limits on the pos- 
sible values of model parameters presented in § [5] We con- 
clude with a summary of our results for the escape fraction 
in § [Uj In our numerical examp les, we adopt the sta ndard 
set of cosmological parameters (iKomatsu et al.l|2009l ). with 
values of f^b = 0.04, = 0.24 and = 0.76 for the 
matter, baryon, and dark energy fractional density respec- 
tively, h — 0.73, for the dimensionless Hubble constant, and 
(78 = 0.81 for the variance of the linear density field within 
regions of radius 8/i~^Mpc. 



2 A GRB DETERMINATION OF THE STAR 
FORMATION RATE DENSITY AT HIGH 
REDSHIFT 

iKistler et al.l (|2009l ) (following lYiiksel et al.ll2008l ) have used 
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the counts of GRBs at high redshift to infer the star forma- 
tion rate density out to z ~ 8.5. Inspection of their Fig- 
ure 4 (see also Figure [T] of this paper) suggests that the 
star-formation rate density is as large as or could even be 
higher at 2: ~ 8 than at z ~ 4 — 6, implying that the ion- 
izing photon emissivity is not falling as observations move 
into the reionization era. This inference requires an assump- 
tion for the initial mass function. However since both ioniz- 
ing photons and GRBs are produced by massive stars, esti- 
mates of the ionizing photon emissivity from the GRB rate 
should be fairly robust against uncertainties in the initial 
mass function (this point is taken up quantitatively in ij |5.5p . 
iKistler et all (|2009l ) argued that we have observed enough 
star formation via high redshift GRBs at z > 6 to reionize 
the universe. This conclusion was reached via a simple esti- 
mate of the number of ionizing photons produced prior to 
z r-u Q given the observed star formation rate density. For 
a Salpeter initial mass function, 4648 ionizing photons are 
produced per baryon incorporated into stars (§[53}. Taking 
this value, together with a constant star formation rate den- 
sity for a time interval At the number of photons produced 
by stars per baryon of IGM is 

^ ^ ( ot) (o.1MqVpc-0 (ioOM^) 

jKistler et al.ll2009l ). The inferred star formation rate den- 
sity at « > 6 implies Af-, ~ 4t2(/osc/0.1). This is con- 
sistent with the modeling of Wyithe & Gen (2007) who 
found Af-y ~ 4 for reionization concluding a.t z — 6 under 
a range of assumptions for the redshift range and efficiency 
of Population-Ill star formation, provided that the escape 
fraction of ionizing radiation is of order 10%. Thus, calcula- 
tions based on the GRB-derived star formation rate density 
imply that there was sufficient star formation to complete 
reionization even with a relatively small value of /esc. 

In this paper we e xtend the simple estimates presented 
m iKistler et al] (HoP). Specifically, we constrain the pos- 
sible values of escape fraction by modeling the reionization 
history and comparing with the observed star formation rate 
density as traced by high redshift GRBs, the ionization rate 
as traced by the Ly-a forest at « < 6, and the optical depth 
to Thomson scattering of CMB photons. The GRB-derived 
star formation rate density values that were presented in 
IKistler et all (j2009l l are special in three ways with respect 
to breaking the degeneracy between escape fraction and star 
formation efficiency during the reionization era. 

Firstly, the star formation rate density is based on num- 
ber counts of GRBs which yields the absolute death/birth 
rate of massive stars. These number counts are independent 
of the escape fraction that can influence other star formation 
rate indicators based on UV emission. The utility of count- 
ing events marking the death of individual stars in order to 
get the star fo r matio n rate density has also been argued by 
iLien fc Fields! (|2009h in the context of supernovae at lower 
redshift. 

Secondly, GRBs are thought to sample the full lumi- 
nosity range of galaxies, so that no correction for luminos- 
ity function is required. In conventional studies of the star 
formation rate, one attempts to correct for both problems , 
but there are uncertainties (e.g.. Hopkins fc Beacomll2006t ). 
The GRB results from .Kistler et al, (2009l ) confirm the two 
large corrections (for obscuration and luminosity function) 



applied to the raw results of iBouwens et al.l l|2008h . It could 
therefore be argued that this agreement implies the com- 
bined corrections are approximately correct, suggesting the 
absence of any deep misunderstandings about the properties 
of high redshift galaxies. Conve rsely, if the lum i nosity func- 
tio n correction applied to the iBouwens et"all (j2008l ) data 
bv IKistler et all (|2009l ') is appropriate, then most of the to- 
tal star formation during the reionization era was within 
small galaxies. With regard to this point it is important to 
note that in more local s amples, most GRB hosts are found 
to be low mass galaxies (jFruchter et al.l I2OO6I : IStanek et al.l 
l2006h . This implies that GRBs are an efficient probe of low 
mass galaxy formation, and is in contrast to studies of the 
UV luminosity function at z > 6 that currently probe only 
the most massive galaxies. Since low mass galaxies seem to 
be dominant during the reionization era we therefore argue 
that GRBs provide the most important constraints on the 
star formation rate density dur ing reionization, and utilize 
the star formation values from IKistler et al.l (|2009| ) in our 
parameter constraints. 

The third fundamental advantage of using GRBs to de- 
termine the star formation rate density is that the data now 
goes deep into the reionization era. This makes them a rather 
direct probe of the star formation responsible for reioniza- 
tion. Importantly, we note that even though the numbers 
of GRBs used in the estimate are low, strong statements 
can still be made about the star for mation rate density 
jKistler et al.ll2009l : lYiiksel et all 120081 ). Of course the ad- 
dition of more high redshift GRBs in the future will aid our 
understanding of the reionization era by lowering the statis- 
tical errors on estimation of the star formation rate density. 



2.1 Summary of star formation rate density from 
high-z GRBs 

As discussed above, IKistler et al.l (|2009l ) have determined 
the star formation rate density deep into the reionization 
era based on the occurrence of GRBs at z > 6. We use as 
constraints the determination of star formatio n rate density 
based on GRBs between 4 < 2 < 8.5 from IKistler et al.l 
(j2009l ). In units of solar masses per year per cubic Mpc, the 
values used are p*,obs ~ 0.061q q5 at z = 7.75±0.75, p^^^obs = 
O.lto.ln at z = 6.5 ± 0.5, p*,obs = O-Olto.ol at z = 5.5 ± 0.5, 
and p*,obs = 0.13^0 05 at z = 4.5 ± 0.5. When coinbined 
with the fact that CMB constraints (|Komatsu et al.ll2009l ) 
place most of the reionization at z ~ 10 this implies that 
we now have multiple constraints probing epochs during the 
reionization era. These data are shown in the upper right 
panel of Figure [T] 



3 ADDITIONAL OBSERVATIONAL 
CONSTRAINTS 

In this section we summarise two additional pieces of data 
that constrain the possible scenarios in addition to the GRB 
derived star formation rate density. 
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Figure 1. Models for the reionization of the IGM and the subsequent post-overlap evolution of the ionizing radiation field. In each 
panel 5 cases are shown, with line-styles corresponding to the cases in Figures 1 2131 The three solid lines shown correspond to the best 
fit (thick line) and 1-sigma deviations for the fiducial case with Tmin = 10-*K and Tion = lO'^K (note that these 3 cases for the fiducial 
model overlap in some of the panels). The dashed line shows a case with Tmin = Tion = 10'*K. The dot-dashed line shows a case with 
^min = lO'^K and Tion = lO^K, but /esc oc (1 -I- z). The models are labeled by the values of /*, /esc and Tes- A value of Ac = 10 was 
chosen for the critical overdensity prior to the overlap epoch. The histories shown correspond to best fit models based on constraints 
from ri2 and pi, a.t z > 4, and also res- Upper Left Panel: The ionization rate as a function of redshift. The observational points are 
fromlBolton fc Haehnel^ bOOTi '). Upper Right Panel: The star formation rate density as a function of redshift. T he data points are from 
iKistler et al Lower Left Panel: The mean-free-path for ionizing photons. The data point is based on Istorrie-Lombardi et al.l 

Lower Right Panel: The volume and mass averag ed fractions of neutral gas in the universe. The observational points for the 
mass- fractions are from the damped Lya measurements of lProchaska et al.l ||2005|) . 



3.1 Thomson scattering optical depth for CMB 
photons 

The measured optical depth to Thomson scattering of CMB 
photons measures the column density of ionized hydrogen 
between the observer and the surface of last scattering. 
The latest observed value of Tes = 0.084 ± 0.016 implies 
that reionization did not occur significantly beyond 2; ~ 10 
l|Komatsu et al.l|2009l ). When computing Tes from the reion- 
ization history we assume that the filling factor of singly ion- 
ized helium equals the ionized fraction of hy drogen, and that 
heliu m becomes doubly ionized at z = 3 (jWyithe fc Loebl 



3.2 Ionization rate from the Lya forest 

Observations of the Gunn-Peterson optical depth in the 
high redshift Ly-a forest imply a very highly ionized IGM 
at 2 < 6, suggesting that reionization was complete by 
that time. The intensity of the ionizing background implied 
by the effective Lya optical depth is sensitive to the de- 
tails of the distribution of gas densities and temperatures. 



but can be reliably mode led via numerical simulation (e.g., 
iBohon fc HaehneltllioOTi '). We take values for the ionization 
rateat 4 < 2; < 6 from the simulations oflBolton fc Haehneld 
(|2007l ). based on the observations of iFan et al. I l|2006l ). In 
units of 10~^^s~^ the values used are Fobs = 0.2'^^^}^ at 
2 = 6, Fobs = 0.5l°:* at 2 = 5, F^bs = ll ag at z =' 4 re- 
spectively. These data are shown in the upper left panel of 
Figure [T] 



4 SEMI-ANALYTIC MODEL FOR 
REIONIZATION 

In this section we summarise the semi-analytic model used 
to calculate the reionization history of the IGM. The basis 
of this model is the differential between ionization and re- 
combination rates for hydro gen in an inhomogeneous IGM. 
iMiralda-Escude et all l|2000l ) presented a model that allows 
the calculation of an effective recombination rate in an in- 
homogeneous universe by assuming a maximum overdensity 
(Ac) penetrated by ionizing photons within HII regions. The 
model assumes that reionization progresses rapidly through 
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islands of lower density prior to the overlap of individual 
cosmological ionized regions. Following the overlap epoch, 
the remaining regions of high density are gradually ionized. 
It is therefore hypothesized that at any time, regions with 
gas below some critical overdensity Ai = pi/ (p) are highly 
ionized while regions of higher density are not. The fraction 
of mass in regions with overdensity below Ai, is found from 
the integral 



FM(Ai)= / dAPv(A)A, 



(2) 



where Pv(A) is the volume weigh t ed pr obability distri- 
bution for A. iMiralda-Escude et al.l l|2000t ) quote a fitting 
function that provides a good fit to the volume weighted 
probability distribution for the baryon density in cosmo- 
logical hydrodynamical simulations. In what follows, we 
draw primarily from the prescription of Miralda-Escude 
et al. (2000) and refer the reader to the original pa- 
per for a d etailed d iscussi on of its motivations and as- 
sumptions. iWvithe fc Loebl l|2003l ) employed this prescrip- 
tion within a semi-anal ytic model of reionization. This 
m odel was extended b y ISrbinovskv fc Wvithd ((20o3) and 
bv lWvithe et al.l l|2008l ). We refer the reader to those papers 
for a full description. 

The quantity Qi is defined to be the volume filling factor 
within which all matter at densities below Ai has been ion- 
ized. The reionization history is quantified by the evolution 
of Qi that evolves according to the rate equation 



dz 



dn-, 



noFmi^i) dz 



aB{l + zf R{/^i)no— + 
dz 



dt dFmiAc) 



dz 



Qi 



i^M(Ai)' 



where as is the case B recombination coefficient, no is the 
comoving density of hydrogen in the mean IGM, and -R(Ai) 
is the effective clumping factor of the IGM. The evolution 
is driven by the rate of emission of ionizing photons per co- 
moving volume dn^/dz. Within this formalism, the epoch 
of overlap is precisely defined as the time when Qi reaches 
unity. Prior to the overlap epoch we must solve for both Qi 
and -Fm (or equivalently Ai). The relative growth of these 
depends on the luminosity function and spatial distribution 
of the sources. In this regime we assume Ai to be constant 
with redshift before the overlap epoch and compute results 
for models with values of Ai = Ac = 10. 

Following overlap Qi = 1 and we describe the post- 
overlap evolution of the IGM by computing the evolution of 
the ionized mass fraction according to the equation 



dfM(Ai) 
dz 



1 dn-j 
no dz 



R(Ai) dt 
Qb 7, — no — . 

a'^ dz 



(3) 



Note that in this post overlap regime the value of Ai is 
the dependent variable describing the ionization state of the 
IGM and is solved for as a function of redshift. 

The emission rate of ionizing photons per co-moving 
volume that is required to compute the reionization history 
can be written 



dz nip dz 



(4) 



where A''.^ is the number of ionizing photons produced per 
baryon incorporated into stars, and is the star forma- 



tion rate per unit volume. As described in the introduc- 
tion, only a fraction of ionizing photons produced by stars 
enter the IGM. Therefore an additional factor of /esc (the 
escape fraction) must be included when computing the ion- 
izing emissivity of galaxies. There are expected to be large 
fluctuations in escape fraction with time and with viewing 
angle for individual gala xies. The escape fra ction could also 
depend on galaxy mass (|Gnedin et al.ll2007l ). However what 
is important for studies of the overall photon budget dur- 
ing reionization is the total cosmic ionizing flux, rather than 
the ionizing fiux near a particular galaxy (although of course 
the latter could affect statistics of Lya absorption or 21cm 
emission from the IGM) . In our model we therefore assume a 
single value for /esc , which should be interpreted as a stellar 
mass averaged value over star forming galaxies. 

The star formation rate per unit volume is com- 
puted based on the coll apsed fraction obtained from the 
iPress fc Schecht"eil (Il974l ') model in halos above the mini- 
mum halo mass for star formation, together with an assumed 
star formation efficiency (/*) 



p* = /^Pc^B (1 - 



m 

di 



rifcol( A/ion) 

dt 



(5) 

where pc is the critical density and the collapsed fraction 
includes separate components from ionized and neutral re- 
gions of IGM. In a cold neutral IGM beyond the redshift of 
reionization, the collapsed fraction should be computed for 
halos of sufficient mass to initiate star formation. The criti- 
cal virial temperature is set by the temperature {Tmin ~ lO** 
K) above which efficient atomic hydrogen cooling promotes 
star formation. Following the reionization of a region, the 
Jeans mass in the he ated IGM limits accretion to halos 
abov e Tion ~ lO'^ K (lEfstathioul 1 19921 : iThoul fc Weinberg! 
|1996| ; iDiikstra et al.|[2004D . Thus, once N-y, Tmin and Tion 
are specified, the reionization model has 2 free parameters 
/* and /esc. 

In order to estimate the ionizing background following 
the end of reionization, our approach is to compute a reion- 
ization history given a particular value of Ac, combined with 
assumed values for the efficiency of star-formation and the 
fraction of ionizing photons that escape from galaxies. With 
this history in place we then compute the evolution of the 
background radiation field due to these same sources. After 
the overlap epoch, ionizing photons will experience attenu- 
ation due to resid ual overdense pockets of HI gas. We use 
the description of iMiralda-Escude et al.l (|2000| ) to estimate 
the ionizing photon mean-free-path [with a reduction of the 
constant of proportionality by a factor of 2 following the dis- 
cussion of Oh fc Furlanctto (2005)], and subsequently derive 
the attenuation of ionizing photons. We then compute the 
flux at the Lyman-limit in the IGM due to sources imme- 
diate to each epoch, in addition to redshifted contributions 
from earlier epochs. 



4.1 Star formation and the production of ionizing 
photons 

In our modeling we assume spectral energy distributions 
(SED) of populati on-II star f orming galaxies using the model 
presented iu iLeitherer et al. I fl999 l. We assume stars of 0.05 
solar metallicity (the effect of metallicity is discussed in 



© 2006 RAS, MNRAS 000, [l]-?? 



6 Wyithe et al. 







o 








Tmln = 10"K Ti„,= 105K - 
(r^^ and f, at z>4) ■ 


CN 




(t^^ ahd.p^ at z>6) 




d 








d 
o 










d 

IN 










O 

d 










o 
d 









0.05 


0.1 

^star 


0.2 




0.05 0. 

^stor 


1 


0.2 


in 








in 








d 








d 












T . =10'^K T 

mm 


= 10*K - 

on 




T . =10'^K T. = 

nnm lon 


lO^K 


f ~(1+z) - 
esc ^ 


CN 

d 






at z>4) - 


c^ 
d 


(r and V, 




z>4) 




c ~ , 














d 








d 








0.05 


: ^^^^ 






m LO 
O 

d 








CN 
O 




^ ^ 




o 








d 








d 








o 
d 








o 
d 










0.05 


0.1 

f , 

star 


0.2 




0.05 0. 

star 


1 


0.2 



Figure 2. Joint constraints on the parameters /* and fesc based on observations of the star formation rate density p*, the ionization 
rate and the optical depth to Thompson scattering of CMB photons Tcs- Upper Left: The fiducial model with Tmin — 

10*K and 

Tion = lO^K. Constraints on and p* are used at z > 4. Also shown for comparison are the constraints on this model if measurements 
of are ignored, which reveals how the GRB star formation data break the degeneracy (gray shading). Upper Right: The fiducial model 
with Tmin = 10*K and Ti^n = lO^K. Constraints on Pi are used at 2: > 6. Lower Left: A model with T^i^ = Ti^n = IC^K. Constraints 
on ri2 and p* are used at 2: > 4. Lower Right: A model with Tmin = 10'*K and Tion = lO^K, and /esc oc (1 + z). Constraints on ri2 
and p* are used at 2; > 4. In each case three contours are shown corresponding to differences in relative to the best fit model of 
= ~ Xmin ~ ^' "^-"^^ ^iid 6.63. Projections of these contours onto the axes provide the 68.3%, 90% and 99% confidence intervals 
on individual parameter values. The large dots show parameters of model histories presented in Figure [T] 



§[S3J, and calculate the number (A^'^,) of ionizing photons 
above the ionization threshold for hydrogen. This number is 
calculated for continuous star-formation, and the low solar 
metallicity is reasonable at high redshift. The mass func- 
tions each have lower and upper limits of O.IM© and 120Mq 
respectively. Taking the functional form dn/dM oc M~" , 
our fiducial model has a Salpeter IMF with a = 2.35 (we 
also discuss other choices for the IMF in § 15.51) . With these 



choices we find A^^, — 4648 in our fiducial model. 



5 CONSTRAINTS ON MODELS OF THE 
REIONIZATION HISTORY 

In this section we confront our model for reionization with 
the three pieces of data described in earlier sections in order 
to constrain the possible scenarios. 

5.1 Parameter constraints 

From the above observations we constrain the possible pa- 
rameterizations of our model for the reionization history. In 
addition to hypothesised parameters Ac, Tmin and Tion, our 



reionization model has two free parameters, /* and /osc, for 
combinations of which we compute the reionization history, 
and calculate the of the model as 



x\f.,U.) = E ^i°sr(/.,/cec)-iogr„.. 



1=0 



or 



E 



logP*(/*,/osc) - logp*,obs 



+ 



'^cs,ob 



(6) 



Here Fobs and p* are the observed ionization rates and star 
formation rate densities measured at a number of redshifts, 
with uncertainty arifi) and a-r (in dex), and Tes^obs is the 
observed optical depth to Thomson scattering of CMB pho- 
tons with uncertainty a-r (in dex). We note that the error 
bars on the observational estimates are not symmetric. 
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5.2 Constraints on the escape fraction in the 
fiducial model 

In the upper left panel of Figure O we present constraints 
on the parameters f-^ and /esc assuming our fiducial model 
with Tmin = lO^K, and Tion = lO^K. As a first illustration 
we have included the available constraints on Tcs as well as 
on r, but have neglected p* in order to highlight the impor- 
tance of its inclusion with the latter examples. The allowed 
region in this case is shown by the gray shading. The ionizing 
background radiation, and the duration of the reionization 
epoch as measured by Te a constrain the combination / star /esc 
tearkan a fc Loebl I2OOII : ISrbinovskv fc Wvithd l2008h . This 
degeneracy can be clearly seen in Figure [2] The dark solid 
contours that are over plotted in the upper left panel of Fig- 
ure |5] show the corresponding constraints following addition 
of measurements of p* from z > 4. The star formation rate 
breaks this degeneracy, since the star formation rate density 
predicted by the model is not directly dependenlQ on /esc- As 
a result, /* and /esc can be separately constrained yielding 
values of 0.085 < /* < 0.145 and 0.035 < /esc < 0.060 for the 
fiducial model. The best fit model has the parameter combi- 
nation (/*, /esc)=(0. 110, 0.046), with a value of Xmin ~ l-l- 
This value is total rather than per degree-of-freedom, and is 
surprisingly low given a 2 parameter fit and 8 data points 
which may indicate that the error bars are overestimated. 
However we note that the values for 3 additional free param- 
eters have been chosen in this model model (Tmin, Tlon and 
Ac). Different values for these are assumed below, leading 
to larger values. 

In Figure [T] we show models (solid lines) cor- 
responding to the best fit, (/*, /csc)=(0. 110, 0.046), as 
well as models near the edge of the 1-sigma contour, 
(/*,/esc)=(0.085,0.060) and (/,, /esc)=(0.145,0.035). The lo- 
cations of these models in the /*-/esc plane are shown by 
dots in Figure [2] In the upper left and upper right pan- 
els of Figure [T] we show the evolution of the ionization rate 
and star formation rate density respectively. In the lower 
left and lower right panels we show the evolution of the 
ionizing photon mean-free-path and mass averaged neutral 
fraction that are also predicted by the model. Note that the 
quantities other than star formation rate density are degen- 
erate for these parameter co mbinations. The data po i nts fo r 
mean- free-path are based oniStorrie -Lombardi et ahl l|l994h . 
The observational points for th e mass-fractions a re fro m the 
damped Lya measurements of iProchaska et al] l|2005h . and 
therefore represent lower limits on the total HI content of 
the IGM. In both cases the curves show excellent agreement 
with these quantities, even though they were not included as 
part of the parameter fit. The values of Tes for these models 
are listed in the lower right panel. The universe is predicted 
to be 50% ionized at z ~ 9.5. 

We note that the observed mean-free-path is found from 
the number density of Ly-limit systems and is independent 
of the Lya forest absorption derived quantities of ionization 
rate and volume averaged neutral fraction, as well as be- 
ing independent of the HI mass-density measurements. Our 
simple model therefore simultaneously reproduces the evo- 

^ There is an indirect dependence of star formation rate density 
on /esc because the star formation rate depends on the fraction 
of universe ionized at a particular time. 



lution of four independent measured quantities, as well as 
the optical depth to electron scattering of CMB photons. 

Despite this success, the accuracy of estimates for the 
ionization rates predicted by the model could be questioned, 
primarily because of the sensitivity of the calculation of the 
mean-free-path to the probability distribution for A. To test 
the sensitivity of our conclusion regarding the value of es- 
cape fraction to the use of F as a constraint, we therefore 
investigate the constraints available when only the obser- 
vations of are used in addition to Tcs- The constraints 
are shown in the upper right panel of Figure (2] In this 
case the best fit model has the same parameter combination 
(/*,/csc) = (0.110,0.046), with a lower value of Xmin = 0.3. 
Thus removing the F points loosens the constraints on /* 
and /esc but does not alter the best fit values. We find val- 
ues of 0.075 < /* < 0.16 and 0.025 < Uc < 0.10 in this 
case. 

In Figure [3] we show the marginalised distribution of 
/esc . The line styles are consistent with those chosen for the 
corresponding contours of in Figure^ In constructing the 
marginalised likelihoods we have assumed prior probability 
distributions for /* and /esc that are flat in the logarithm, 

i.e., dPprior/dlog/* OC 1 and dPpHor/dlog /esc oc 1. 



5.3 Effect of radiative feedback on low mass 
galaxies 

There is some theoretical uncertainty surrounding the 
strength of radiation feedback on the formation of low mass 
galaxies (Dijkstra ct al. 2004; Mc singcr fc Diikstra 2008). 
We therefore repeat our constraints on /* and /esc using 
a modification of the fiducial model with no radiative feed- 
back, i.e., Tinin — Tion ~ lO'^K. The constraints based on 
this model are shown in the lower left panel of Figure (2] 
The best fit modified model has the parameter combination 
(./■*,, fosc)=(0.070,0. 055), with a value of Xmin = 3.0 (this 
model is therefore disfavoured by the data relative to the 
fiducial case). As reionization progresses, more gas is al- 
lowed to accrete into galaxies in this case than in the fiducial 
model. As a result the star formation efficiency required to 
produce the observed star formation rate density is reduced. 
We find values of 0.05 < /* < 0.095 and 0.04 < /esc < 0.08 
for this model. Thus the details of the stellar mass accre- 
tion rate influence the inferred star formation efficiency, but 
only weakly effect the inferred escape fraction of ionizing 
photons. 

Figure [3] shows the corresponding marginalised distri- 
bution for escape fraction (dashed line). As noted above, 
the removal of radiative feedback from the model has little 
effect on the derived escape fraction distribution. On the 
other hand the reionization history is quite different in this 
case. The dashed lines in Figure [T] show the best flt model in 
this case (/*, /osc)=(0. 07, 0.055). The lower value of /* means 
that reionization gets underway later. However despite this 
reionization is finished earlier, with an associated rise in the 
ionizing background at higher redshifts in this case than in 
the fiducial model. This rapid rise occurs because reioniza- 
tion is assumed not to be self limiting in this case and so is 
completed by low mass galaxies. 
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Figure 3. Constraints on the parameter /eac corresponding to 
marginalisation of tiie joint distributions in Figure [2] (the line 
styles are the same for the joint and marginahsed distributions in 
these figures). Solid line: The fiducial model with Tmin = fO^K 
and Tioji = fO^K. Constraints on and p^, are used at 2 > 4. 
Dotted line: Constraints on the fiducial model without inclusion 
of ri2. Dashed Line: A model with T^\^ = Tjon = fO'^K. Con- 
straints on and p* are used at z > 4. Dot-Dashed Line: A 
model with Tmin = lO^'K and Tion = lO'^K, and /esc oc (1 + z). 
Constraints on ri2 and are used at z > 4. 

5.4 Evolving escape fraction 

Our modeling so far has assumed an escape fraction that is 
constant with redshift. To test the impact of this assumption 
on the value of /esc we therefore consider a case where our 
fiducial model is modified so that /csc(z) = /csc(l + z)/7, 
and constrain the value of /esc a.t z = 6. The results are 
shown in the lower right panel of Figure [5] In this case 
we find values which are very similar to the fiducial case, 
0.085 < /* < 0.145 and 0.035 < /esc < 0.06 for this 
model. The best fit model has the parameter combination 
(/*, /csc)=(0. 11, 0.053), with a value of Xmin ~ 3-1 (this 
model is therefore also disfavoured relative to the fiducial 
case) . 

Figure [3] shows the corresponding marginalised distri- 
bution for escape fraction (dot-dashed line). The evolution 
of escape fraction is degenerate with evolution in the mass 
accretion rate of galaxies, and so as in the previous example 
has little effect on the derived escape fraction distribution. 
The dot-dashed lines in Figure [T] show the best fit model 
in this case (/*, /csc) = (0. 11, 0.053). The reionization history 
is again quite different from the fiducial model. The larger 
value of /esc towards high redshift results in an earlier reion- 
ization and a higher Tcs- 

5.5 Uncertainty in the choice of IMF 

Both the number of ionizing photons produced per star as 
well as the star formation rate density inferred from the 
number of GRBs observed during the reionization era de- 
pend on the IMF assumed. The former dependence arises 
because more massive stars are hotter and therefore emit 
a larger fraction of their energy at frequencies above the 
ionization threshold for hydrogen. Regarding the latter de- 
pendence, we note that the relation between the GRB rate 
and the star formation rate density is calibrated at redshifts 
below 2 ~ 4 where they can each be separately determined 



iKistler et al.ll2009t ). Most of the star formation rate indica- 
tors (e.g., Ha) probe only massive stars, and hence a cor- 
rection involving an IMF is required. However, what is im- 
portant for the present analysis is the number of ionizing 
photons implied by the observed number of high redshift 
GRBs, i.e., relating one quantity depending just on massive 
stars to another. With respect to this quantity the depen- 
dencies on IMF mentioned tend to cancel, because a more 
top-heavy IMF leads to increases in both the GRB rate and 
the number of ionizing photons produced. We note here that 
within the context of studies of the star formation history, 
top-heavy can really be interpreted to mean bottom-light. 
This is because one measures the light (or the birth/death 
rate) from the massive stars and infers the contribution to 
stellar mass from lower masses. In other words the analysis 
is normalized by the high-mass end of the IMF. This situa- 
tion would be reversed if one were looking at the integrated 
stellar mass, which is dominated by low-mass stars. 

To assess the systematic error in determination of the 
escape fraction introduced by the choice of IMF we consider 
a range of initial mass functions (IMF), as summarised in 
Table [1] As in our fiducial case, we assume stars of 0.05 so- 
lar metallicity, and present the number of ionizing photons 
above the ionization threshold for hydrogen. This number 
is again calculated with mass functions that have lower and 
upper limits of O.IM© and I2OM0 respectively. We consider 
two cases with turnover in the mass function below O.5M0 
(with an index in that mass range of a = 1.5). These cases 
have high mass in dexes of a = 2.35 (lab eled Salpeter A) 
and a = 2.15 (jBaldrv fc Glazebrookl l2()03i . labeled BG) re- 
spectively. The Salpeter A and BG IMFs were identified by 
iHopk ins fc BeacomI (I2OO6I ) as providing reasonable bounds 
on the normalisation of the cosmic star formation history. 
We also consider the case of a top-heavy mass function hav- 
ing a single index of q = 1.95 across the entire mass range. 
The values for A'^, are listed in Table [T] demonstrating that 
more top-heavy mass functions produce a larger fraction of 
their energy as photons above the ionization threshold. The 
ratio relative to the value A'sai = 4648 for a Salpeter mass 
function is also presented in the second last column for ease 
of comparison. 

In the 4th column of Table [T] we show the fraction of 
the star formation inferred assuming a Salpete r IMF that 
would have been derived in iKistler et all l|2009l ) if a differ- 
ent IMF had been used instead. Where a more top-heavy 
IMF is assumed, the star-formation rate density implied by 
the observed GRBs is reduced by this factor relative to our 
fiducial model. Since a change in the assumed IMF leads to 
a change in the inferred star formation rate density, it also 
leads to a proportional change in the derived star formation 
efficiency. Thus for each IMF listed the value of p*/ Pazi rep- 
resents the size of the systematic uncertainty in constraints 
on /*. 

The change in IMF also leads to a change in N-y, and 
the escape fraction required to reproduce the observables of 
reionization is therefore proportional to {f)i,N^)~^ . In the 
final column we list the quantity 

/iMF = P*/pssl X N^/N^^i, (7) 

which gives the ratio of ionizing photons to GRBs, relative 
to the Salpeter IMF. The range in values of /imp provides an 
estimate of the systematic uncertainty in the derived escape 
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Table 1. Summary of IMFs used in this work. 



IMF 


"low 


a 


P*/Psal 


N-y 




/imp 


Salpeter 




2.35 


1 


4648 


1 


1.00 


Salpeter A 


1.5 


2.35 


0.77 


6034 


1.30 


1.00 


BG03 


1.5 


2.15 


0.50 


11139 


2.39 


1.20 


Top-Heavy 




1.95 


0.37 


17553 


2.75 


1.39 



> 10%, ionizatio n by secondary electrons becomes inefficient 
l|Ricotti fc Ostriker 2004). As a result, X-ray ionization is 
self regulating at the level of 10 — 20%. Modeling the 
possible contribution of X-rays to reionization is beyond the 
scope of this paper. However the self-regulation implies that 
neglect of X-rays would lead to an overestimation of /esc 
by this factor at most, and we therefore estimate a system- 
atic component of uncertainty owing to the unknown X-ray 
contribution to reionization of ~ 0.05 dex. 



fraction owing to the choice of IMF. As expected, the net 
effect on the reionization considerations caused by choosing 
a different IMF is quite modest. It is possible that the IMF 
is evolving with look back time. In this case, by extension 
of the above arguments our estimate of the escape fraction 
would sustain only a small systematic uncertainty. Impor- 
tantly, we find that for metal enriched stellar populations, 
the choice of IMF contributes a systematic uncertainty on 
the derived escape fraction of only ~ 0.1 dex. 

Finally, we note that the assumed metallicity can have 
an effect on both the calibration of star formation rate den- 
sity, and on the value of A''.^ computed for a particular IMF. 
Calibrations of star formation rate density are performed 
locally and so refer to solar metallicity. On the other hand 
we have assumed a metallicity of l/20th solar in comput- 
ing the reionization history, as may be more appropriate at 
high redshift. The Salpeter IMF with solar metallicity yields 
N.y = 3328 ionizing photons per baryon, a factor of 1.4 times 
smaller than in our fiducial model. This metallicity uncer- 
tainty corresponds to an additional ~ 0.1 dex uncertainty 
on /esc- It may be possible to remove this metallicity un- 
certainty in the determination of /esc by noting that many 
star formation rate indicators (e.g.. Ha) are proportional to 
the ionizing flux, so that the higher number of ionizing pho- 
tons for low metallicity stars is offset by a smaller inferred 
star-formation rate. We have not pursued this in the cur- 
rent work. We conservatively estimate that accounting for 
the uncertainties in IMF and metallicity leads to ~ 0.2 dex 
of systematic uncertainty in the determination of /esc from 
our analysis. 



5.6 Reionization by X-rays 

Our analysis has thus far assumed that reionization is en- 
tirely due to high mass stars. However in addition to UV 
photons from the first galaxies and quasars, it has been 
suggested that a background of X-ra y photons at very 
high redshift may be important (e.g., iRicotti fc Ostrikeil 
I2OO4I '). Indeed, several authors have proposed an initial 
phase of preheating and partial ionization of the IGM 
by X -rays, resulting primarily fr o m black hole accretion 



Shull fc van Steenber 3 Il985l : IVenkatesan et all I2001I : 



iRicotti fc Ostriker! I2OO4I ). An X-ray background could also 
have originated from X-ray binaries and supernova rem- 
nants. In difference to UV photons. X-rays ionize hydrogen 
both directly and through secondary ionizations by photo- 
electrons from ionized helium, with the latter dominating 
so that many hydrogen atoms can be ionized by a single 
X-ray photon. As reionization proceeds, however, the pro- 
portion of each X-rays energy that is deposited into the 
IGM as heat, increases. In particular, for ionization fractions 



5.7 Numerical constraints on the escape fraction 

Our fiducial model constrains the escape fraction to have a 
value log /esc = — 1.35±0.15 (68%). Of the systematic errors 
discussed the largest arise due to uncertainty in the IMF and 
in the metallicity of the star forming populations responsible 
for reionization. We estimate that each of these contribute 
^ 0.1 dex of systematic uncertainty. Thus in summary of 
our analysis we find log /esc ~ ( — 1.35 ± 0.15) ± 0.2. 



6 SUMMARY 

The fraction of ionizing photons that escape their host 
galaxy and so are able to ionize hydrogen in the inter- 
galactic medium (IGM) is a critical parameter in studies 
of the reionization era and early galaxy formation. Indeed 
as can be seen directly from equation ([l}, the amount of 
star formation required to reionize the Universe is inversely 
proportional to the escape fraction. Determination of the 
expected value for the escape fraction is problematic for in- 
dividual galaxies. Observationally, very deep spectra must 
be obtained to overcome absorption in the IGM and detect 
photons that escape the host galaxy at energies beyond the 
Ly limit. Theoretical calculation of escape fraction is diffi- 
cult owing to the large dynamic range and complex gas and 
star formation physics that must be modeled to resolve the 
clumpy ISM. 

Many studies of the reionization history have compared 
theoretical models to the electron scattering optical depth 
of CMB photons, and the ionization rate from the Lya for- 
est at the end of the overlap era. These studies are sub- 
ject to degeneracies between the escape fraction of ioniz- 
ing photons, and the efficiency wit h which high redshift 
galaxies turn mass into stars (e.g. , IWvithe fc Loebl l2003l: 
llaiman fc Holder 2003 : ICenI l2003l : Ichoudhu rv fc Ferraral 
I2OO5I : IShull fc Venkatesann2008l '). While upper limits have 
previously been obtained for surveys of high redshift star 
forming galaxies (|Bouwens et al.|[2008l ). the recent determi- 
nation of the star formation rate density during the reion- 
ization era from the most distant Gamma Ray Bursts dis- 
covered provides the first method for directly measuring the 
star formation rate density during the reionization era. With 
respect to calculating the escape fraction, this method has 
the advantages over estimates of star formation rate from 
the galaxy luminosity function that it directly observes the 
total star formation rate (i.e., no correction for a galaxy flux 
limit), and is well calibrated at lower redshift (jKistler et al.l 
i2009). 

In this paper we have used this first determination of 
the star formation rate density at 2: > 6 to break the degen- 
eracy between star formation efficiency and escape fraction. 
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Our analysis employed a semi-analytic model for reioniza- 
tion that we confronted with three complementary sets of 
observations, i) the electron scattering optical depth of CMB 
photons, ii) the ionization rate from the Lya forest at the 
end of the overlap era and Hi) The star formation rate den- 
sity out to z = 8.5. By constraining the model parameters 
with these observations we show that the escape fraction of 
ionizing photons from high redshift galaxies is /esc ~ 5%. 
For our fiducial model we find log /esc = —1.35 ± 0.15 

Importantly, this value is stable against a range of sys- 
tematic uncertainties in the modeling, including the imple- 
mentation of radiative feedback and the choice of IMF. In 
agreement with previous studies we find that the mass as- 
sembly history of galaxies was important for the reionization 
history, and impacts on the inferred value of star formation 
efficiency. However, the value of escape fraction inferred is 
not sensitive to the details of the mass assembly. We also 
investigated an evolving escape fraction, finding that this 
did not strongly influence the value of escape fraction in- 
ferred for galaxies at z = 6. Our modelling does make the 
assumption that UV photons produced by stars dominate 
the reionization process. A possible cont ribution of X-rays 
to reionization (jRicotti fc Ostrikerl |2003 ) complicates this 
interpretation, and implies that our escape fraction could 
be overestimated by a factor up to 1.1 — 1.2. Both the num- 
ber of ionizing photons and the number of GRBs per unit 
star formation are sensitive to the IMF. The dependencies 
on IMF tend to cancel in our analysis, however we consid- 
ered a range of possibilities. We estimate an uncertainty of 
^0.1 dex in /esc owing to the unknown shape of the IMF, 
and another ~ 0.1 dex owing to the unknown metallicity 
of the galaxies responsible for reionization. Including these 
uncertainties we express our estimate of escape fraction as 
log /esc = ( — 1.35 ± 0.15) ± 0.2. Our analysis assumed that 
the escape fraction is independent of galaxy mass. Previous 
work has sugge sted that lower mass galaxies have smaller es- 
cape fractions (|Gnedin et aLlbOOTh . In this scenario our es- 
cape fraction should be considered as a stellar mass weighted 
average over star forming galaxies. 

Our simple semi-analytic model is able to provide a 
statistically acceptable fit to the available observations un- 
der a range of assumptions. This implies that models of 
the reionization history are under constrained [although 
more variables could potentially be used as constraints 
IChoudhurv fc Ferraral (|2005l . l2006l ')]. In the near future 21cm 
observations of neutral IGM d uring the r e ioniza tion era will 
provide a powerful new probe. iLidz et al] l|2008l ) have shown 
that the first generation of low frequency telescopes will 
be able to measure the neutral fraction accurately during 
epochs when the universe is ~ 50% ionized. The exam- 
ple reionization histories presented in this and other pa- 
pers illustrate that this will provide important discriminates 
among possible reionization histories. With respect to the es- 
cape fraction, such measurements in combination with the 
evolution in star formation rate density will allow for the 
evolution in addition to the value to be measured. 
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